Dynamical cluster approximation study of the anisotropic two-orbital Hubbard model 
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We investigate the properties of a two-orbital Hubbard model with unequal bandwidths on the 
square lattice in the framework of the dynamical cluster approximation (DCA) combined with a 
continuous-time quantum Monte Carlo (CT QMC) algorithm. We explore the effect of short-range 
spatial fluctuations on the nature of the metal-insulator transition and the possible occurrence of 
an orbital-selective Mott transition (OSMT) as a function of cluster size N c . We observe that for 
N c = 2 no OSMT is present, instead a band insulator state for both orbitals is stabilized at low 
temperatures due to the appearance of an artificial local ordered state. For N c = 4 the DCA 
calculations suggest the presence of five different phases which originate out of the cooperation and 
competition between spatial fluctuations and orbitals of different bandwidths and an OSMT phase 
is stabilized. Based on our results, we discuss the nature of the gap opening. 

PACS numbers: 71.10.Fd,71.27.+a,71.30.+h,71.10.Hf 



The correlation driven metal-insulator transition 
(MIT) in two-dimensional (2D) correlated systems is still 
poorly understood. While the behavior of 2D one-band 
systems at half filling seems to be settled, this is not 
the case for multi-orbital systems. In one-band systems 
long range correlations or local order in small size clus- 
ters caused by perfect nesting (Slater physics) [E 0, 0, 0] 
open a gap in the weak-coupling regime, whereas in the 
strong-coupling regime the on-site Coulomb repulsion is 
the driving force for the gap opening (Mott physics) Q. 
On the other hand, even though multi-orbital models are 
better suited to describe real materials @ and display 
rich phase diagrams @,@,|{| as well as interesting phy sics 
like the orbital-selective Mott transition (OSMT) pJOlfTll ]. 
they are still under debate. The reason for that is their 
more complex structure compared to the one-band model 
due the orbital degrees of freedom, crystal field splitting- 
effects and the Hund's rule exchange coupling. This sit- 
uation gives us a strong motivation to investigate the 
two-orbital system. 

The OSMT has recently been intensively studied in 
the context of a weakly correlated band coexisting and 
interacting with a more strongly correlated one in a two- 
orbital system [13, [H, E3, E EE E3, H, El, ■ Is- 
sues like (1) the importance of full or Ising-type Hund's 



rule coupling [12], [lj, |16|, [18|, (2) the consequences of 
anisotropic Hund's rule coupli ng 11911 . (3) the role of the 
ratio of the two bandwidths [la], (4) the inclusion of 
the hybridization between bands [lj], (5) the effect of 



crystal field splitting [17[ and (6) the extension to the 
three-band case [2l| have already been addressed. Nev- 



ertheless, the importance of spatial fluctuations has not 
yet been explored since most calculations have been per- 
formed within the single-site dynamical mean field the- 
ory (DMFT) 13, 23| where spatial fluctuations are com- 



pletely ignored. On the other hand, it has already been 
noticed that even in a single-band case, inclusion of spa- 
tial correlation will qualitatively change the scenario of 
the Mott MIT @, E3, El] ■ Therefore, it is crucial to ad- 
dress the effect of spatial fluctuations on the OSMT and 
the phase diagram. 

Very recently, Bouadim et al. [H studied the OSMT 
by means of a determinant quantum Monte Carlo method 
(DQMC) on the square lattice and showed that an itin- 
erant band can coexist with a fully localized band in a 
two-orbital Hubbard model as long as long range an- 
tiferromagnetic correlation is absent. However, since 
the DQMC calculation was based on a simplified model 
where one of the two orbitals is constrained to be fully 
localized, it still remains unclear whether the OSMT 
survives in the system with spatial fluctuations or not. 



Moreover, since previous DMFT [H, [iH 111 fl5L [T_ 
18. 19] and a slave spin mean field calculation 21] arc 
based on the Bethe lattice, it is interesting to move in the 
direction of real systems by studying the case of a two- 
dimensional model on the square lattice with the Fermi 
level at a van Hove singularity at half filling. 

In this Letter we concentrate on the nature of the 
gap opening and the OSMT in a two-dimensional sys- 
tem. The anisotropic two-orbital Hubbard model on the 
square lattice at half-filling is the simplest model which 
can describe the OSMT including spatial fluctuations. 
The Hamiltonian is given as 

H ^ ^ lmCi 7niJ Cjm<T "F U ^ ^ ^tm|^™| 

(ij)ma irn 

+ V ,(U' - Jz)naani2<7>, 

* ' 1(7(7 

where t m for orbital m = (1, 2) denotes the hopping in- 
tegrals between nearest-neighbor (n.n.) sites i and j, U 



2 



and V are intra-orbital and inter-orbital Coulomb repul- 
sion integrals respectively and J z nnani2a for spin a is the 
Ising-type Hund's rule coupling term. In our calculations 
we ignore spin-flip and pair-hopping processes. We also 
set ti/t — 0.5 (narrow band), t 2 /t = 1.0 (wide band), 
J z = {7/4 and U' = U/2. For this model we employ the 
dynamical cluster approximation (DCA) method with 
cluster sizes N c = 2 and 4. The DCA method [USSS! 
can not only overcome the problem of the single-site 
DMFT method [H SI , where Mott physics rather than 
Slater physics is emphasized in the paramagnetic phase 
due to the lack of spatial fluctuations, but it is also com- 
putationally cheaper than lattice calculations. We use a 
weak coupling continuous-time quantum Monte Carlo al- 
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gorithm as an impurity solver [2jJ,|30j. We shall present 
results on the spin-spin correlations, double occupancy, 
self-energy and density of states (DOS). 

First, let us discuss the results obtained from the DCA 
with a two-site cluster (N c = 2). It is known that, for a 
one-band system at low temperatures, the formation of a 
local singlet state driven by Slater physics is responsible 
for the gap opening. In the two-band system orbital fluc- 
tuations are present. Due to the Hund's rule coupling and 
the Coulomb interaction, ferromagnetic (FM) correlation 
between orbitals and antifcrromagnetic (AF) correlation 
between sites develop. In order to check for these cor- 
relations we measure the on-site (inter-site) inter-orbital 
spin-spin correlations (Sf-^S^) ({^ii^i+i 2)) as a func- 
tion oi U/t. The results are shown in Fig. [Ha) . 

As the Hund's rule coupling J z = U /4 is increased, the 
on-site inter-orbital FM correlations (positive sign) and 
inter-site inter-orbital AF correlations (negative sign) are 
enhanced for both temperatures T/t = 0.3andT/t = 0.1. 
At high temperatures T/t = 0.3 the on-site inter-orbital 
FM correlations are stronger than the inter-site inter- 
orbital AF correlations which can be attributed to the 
fact that thermal fluctuations suppress the AF correla- 
tions. At T/t = 0.1 both correlations are of the same 
magnitude. This behavior suggests the appearance of a 
local ordered state in the low temperature regime. In 
order to verify whether this state, in analogy to the one- 
band model, is responsible for the gap opening as de- 
scribed by Slater physics, we calculate the temperature 
dependence of double occupancy for both orbitals. If 
Slater physics is dominant, as the temperature is de- 
creased the formation of the local order which reduces 
the potential energy U {n-\n^) should cause the gap open- 
ing. In Fig. QJb) we present the double occupancy as a 
function of temperature T/t for U/t = 2.4. The double 
occupancy in both bands decreases with decreasing tem- 
perature and it shows a more abrupt drop near T/t = 0.2. 
This behavior gives strong evidence of Slater physics, and 
the band insulator in both orbitals should be present 
at zero temperature for all positive interaction strengths 
U/t. 

Next, we explore the DCA for a 4-site cluster (N c = 4). 
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FIG. 1: (a) On-site and inter-site inter-orbital spin-spin cor- 
relations for N c = 2 as a function of U/t at temperatures 
T/t = 0.1 and T/t = 0.3. (b) Double occupancy for N c = 2 
as a function of T/t for U/t = 2.4. Left and right axes are for 
the double occupancy of the narrow band and the wide band, 
respectively. 



The inclusion of next nearest-neighbor (n.n.n.) correla- 
tions in N c = 4 suppresses the local ordered state en- 
hanced artificially in the N c = 2 cluster. In addition, the 
system shows a weak degree of frustration because of the 
absence of long range correlations. In this way, the Mott 
physics present in the single-site systems coexists with 
Slater physics present in the two-site systems. Therefore 
we believe that the description in terms of the N c — 4 
clusters is closer to the real materials at finite tempera- 
ture. In Fig. HJa) we compare the on-site and inter-site 
inter-orbital spin-spin correlation results for N c = 2 and 
4 for T/t = 0.1. The same magnitude of both corre- 
lations for N c = 2 implies the presence of a relatively 
strong inter-site local ordered state, while the deviation 
of those for N c = 4 indicates that the local ordered state 
is released due to the influence of the n.n.n. correla- 
tions. In order to compare directly the strength of this 
local ordered state for N c = 2 and 4, we plot n.n. and 
n.n.n. correlations for the narrow band in Fig. EJb). In 
the weak coupling regime the n.n. correlation strength 
for N c — 2 due to enhanced Slater physics is larger than 
that for N c = 4. In the strong coupling regime the n.n. 
correlation strength is similar for N c = 2 and 4 since the 
insulating state for N c = 4 is induced by cooperation 
of Mott and Slater physics. We also find strong n.n.n. 
correlations for N c = 4. 

The competition among magnetic and orbital fluctu- 
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FIG. 2: (a) The on-site (inter-site) inter-orbital spin-spin 
correlations for N c — 2 and 4 as a function of f//t at 
T/t — 0.1. (b) The narrow band (next) nearest-neighbor 
spin-spin (n.n.n.,n.n.) correlations for N c — 2 and 4 as a 
function of U/t at T/t = 0.1 



ations as well as weak frustration for N c = 4 should 
generate a rich phase diagram. In order to investigate 
this complex situation we analyze, in what follows, the 
on-site self-energy. The imaginary part of the on-site self- 
energy ImS(iw n ) provides information about the possi- 
ble Fermi- liquid/non- Fermi- liquid behavior of the system 
as well as the nature of the gap opening. In Figs. [3Ja) 
and (b) we present ImE(ia; n ) for the narrow and wide 
band, respectively, at T/t = 0.1. According to Fermi- 
liquid theory, ImS(w) at T = at to — > extrapolates 
to 0. In the weak-coupling regions below U/t — 1.4 this 
Fermi-liquid behavior is seen in both bands. Between 
U/t = 1.4 and 1.8 Fermi-liquid behavior is still present 
in the wide band, while non-Fermi-liquid behavior is ob- 
served in the narrow band. The electrons begin to local- 
ize in the narrow band driven by both Slater and Mott 
physics, while those in the wide band are still delocalized. 
As the interaction is increased, non-Fermi-liquid behav- 
ior is observed in both bands. At U/t = 2.8, Im£(iw n ) 
in the narrow band diverges, which indicates the open- 
ing of a gap, while the metallic state (non-Fermi-liquid) 
is still present in the wide band. These results evidence 
a OSMT. In the strong-coupling region U/t = 3.4 the 
insulating state is observed in both bands. 

In what follows we shall analyze the nature of the 
gap opening. According to recent results obtained for 
the single-band plaquette Hubbard model [HI, momen- 
tum sectors K=(Q, 0)/(tt, 7r) and (w, 0)/(0, 7r) undergo 
a metal to band insulator transition and a metal to 
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FIG. 3: The imaginary part of the on-site self-energy for 
U/t = 1.0,1.4,2.2,3.0,3.4 at T/t = 0.1 for N c = 4 (a) in the 
narrow band and (b) in the wide band. The real (c) and 
imaginary (d) part of self-energy at the lowest Matsubara fre- 
quency Wo for K=(n,n) and (n,0) sectors as a function of 
U/t. 



Mott insulator transition, respectively. In Figs. [3] (c) 
and (d), respectively, we present the real and imagi- 
nary parts of the self-energy at the lowest Matsubara 
frequency cl>o, ReE(kjn) and ImS(iwo)) for K=(tt 1 tt) 
and (7r,0) in both bands. While ReE(iwn) gives infor- 
mation about the energy shift of the spectral function, 
ImE(iwo) introduces the scattering rate. As the interac- 
tion is increased, ReS^^) = — Re£( ,o) increases while 
ImE( W)7r ) = ImS( 0) o) remains small in both bands. These 
results suggest a metal to band insulator transition where 
the gap is opened through separation of the poles away 
from the Fermi level. On the other hand, as the interac- 
tion increases ImE^ o) = ImS( 0i7r ) displays a divergent 
behaviour and ReE( ff- o) = ~~ ReE(n i7I -) in both bands is 
zero due to particle-hole symmetry. Therefore, in the 
strong coupling region, the gap in the K—(tt, 0) and (0, w) 
sectors is only induced by the divergence of ImE(ia;n) 
which is a signature for Mott physics. These results 
are similar to the single-band plaquette Hubbard model 
results [25l | but, while a first-order transition occurs in 
the single-band Hubbard model, the OSMT behaviour 
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FIG. 4: Density of states at T/t = 0.1 and N c = 4 for (a) 
U/t — 2.8 and (b) U/t = 3.4. We employ the Pade approxi- 
mation method for the analytic continuation. 
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FIG. 5: The phase diagram with five phases for N c — 4. 



is present in the two-band Hubbard model. In order to 
show the OSMT more clearly, we present in Fig. 0] the 
density of states (DOS) at T/t = 0.1 for the interaction 
values U/t = 2.8 (Fig. g] (a)) and U/t = 3.4 (Fig. g] (b)). 
These interaction values have been identified as onsets 
for the OSMT and insulator phases, respectively. For 
U/t = 2.8 the narrow band exhibits a gap at the Fermi 
energy (ui = 0), while the wide band has a finite DOS at 
u> = 0. This means that for a given interaction strength 
U/t two stages of the Mott transition are present, with 
a Mott insulator in the narrow band and a metal in the 
wide band. At U/t = 3.4 both bands show a gap at 
lu = 0. The gap in the narrow band is wider than that in 
the wide band. Finally, we plot the phase diagram with 
the identified five phases for N c — 4 in Fig. El 

In summary, we have explored the anisotropic two- 
orbital Hubbard model using the DCA method with clus- 
ter sizes N c = 2 and 4. The DCA cluster with N c = 2 
for the single-band model is known to describe a system 
with artificially strong local order between sites and the 
gap opening is controlled by Slater physics. Our results 
show that this inter-site AF correlation is still strong in 
spite of orbital fluctuations, leading to a gap at low tem- 
peratures. The appearance of the insulating states can 
be described by Slater physics. We have also investi- 
gated within DCA the N c — 4 cluster which includes 
n.n.n. correlations. Unlike the N c = 2 cluster, the lo- 
cal ordered states are not present in the weak-coupling 
limit. In the very weak-coupling regime Fermi-liquid be- 
havior is present in both bands. As the interaction in- 



creases, the electrons in the narrow band weakly localize 
and non-Fermi-liquid behavior is observed, even though 
the Fermi-liquid behavior is still present in the wide band. 
In the intermediate region, non-Fermi-liquid behavior is 
observed in both bands. In the strong-coupling region the 
electrons in the narrow band are completely localized and 
those in the wide band are partially localized which can 
be described as the OSMT. In the very strong-coupling 
region both orbitals are insulating. The nature of the 
gap opening is that of coexisting Slater physics in the 
momentum sector K=(0,0)/(tt,tt) and Mott physics in 
momentum sector K=(ir, 0) / (0, it). 
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